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Abstract 

Structure of eight superdeformed bands in the nucleus ^^^Tb is analyzed using 
the results of the Hartree-Fock and Woods-Saxon cranking approaches. It is 
demonstrated that far going similarities between the two approaches exist and 
predictions related to the structure of rotational bands calculated within the 
two models are nearly parallel. An interpretation scenario for the structure 
of the superdeformed bands is presented and predictions related to the exit 
spins are made. Small but systematic discrepancies between experiment and 
theory, analyzed in terms of the dynamical moments, J'^'^\ are shown to exist. 
The pairing correlations taken into account by using the particle-number- 
projection technique are shown to increase the disagreement. Sources of these 
systematic discrepancies are discussed ~ they are most likely related to the 
yet not optimal parametrization of the nuclear interactions used. 
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I. INTRODUCTION 



In the present article the results of microscopic calculations for the high spin properties of 
the superdeformed (SD) rotational bands in the nucleus ^^^Tb are presented. This nucleus is 
a direct neighbor of the "doubly magic" superdeformed nucleus ^^^Dy which is characterized 
by particularly large gaps in the neutron and proton deformed single-particle spectra. We 
investigate the rotational bands together with their microscopic structure in terms of the 
single-particle orbitals. 

The theoretical description of such properties, most often applied in the literature, is 
based on the average field approximation (see Ref. for a recent review). Most of the 
present knowledge about the nucleonic structure underlying the properties of the SD rota- 
tional bands comes from analyzing the single-particle Routhians generated by making use 
of the rotating Nilsson or deformed Woods-Saxon (WS) Hamiltonians. The total nuclear 
energy as a function of spin is then obtained by using the Strutinsky total energy expression 
(see for example Ref. and references therein). Similar analyses based on the Hartree-Fock 



(HF) 0-0, or relativistic mean field (RMF) p!0|-[T^ approximations are more difficult and 
have not been, up to now, applied for global calculations, with the notable exception of Ref. 



I^ . However, they have an advantage of directly providing the total energy as a function 
of spin which can be compared with the experimental data. In the present work we apply 
and compare the WS and HF cranking approximations. 

Uncertainties in predictions coming from various models lead to ambiguities in the in- 
terpretation of the experimental results. The problem of a (non-) uniqueness in the inter- 
pretation of the microscopic structure of SD bands is strongly aggravated by the fact that 
the spins of nearly all measured bands in question are not known. Moreover, the nuclear 
interactions responsible for the single-particle level ordering in the strongly deformed nuclei 
are obtained mainly by extrapolating the properties of the interactions known from spherical 
and moderately deformed nuclei. 

Direct tests of the structure of wave functions of the single-nucleonic orbitals (via e.g. 
magnetic moments) are considered extremely difficult from the experimental point of view 
due to the very short life-times of the SD high-spin states. In such a situation we are 
confronted not only with the problem of a microscopic interpretation of the quickly growing 
body of the experimental information — but also with the problem of a cross-check analysis 
of theoretical predictions coming from distinct nuclear interactions used. It will be one of 
our objectives to provide such a test by comparing the microscopic properties calculated 
within the Strutinsky and HF methods. 

For a realization of our objective, the experimental knowledge of eight bands in one single 
SD nucleus, ^^^Tb, is extremely useful. At present there are only very few nuclei in which 
the number of known SD bands is comparable to that in ^^^Tb. Moreover, the transition 
energies of those bands behave regularly as functions of the angular momentum /. Such 
a regularity provides the optimal conditions for theoretical verifications of single-nucleonic 
configurations via the dynamical moments defined by 

J^^ = ^. (1.1) 

auj 

This is so, because the abrupt band crossings (relatively abundant in other nuclei) may 
otherwise cause strong variations of jT'-^^-moments as a function of cu, thus rendering the 
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comparisons of the typical behavior (and values) of the JT*^^) -moments with experiment less 
direct. In some cases the band crossings may play an important role in helping to identify 
some orbitals via the characteristic crossing frequencies. However, we found it preferable to 
first discuss the structural elements, such as the global behavior of the jT^^^-'-moments and 
their values, in relation with some characteristic nucleonic configurations. Another group 
of problems related to the band crossings, crossing frequency systematics and alignments is 
left for future investigations. 

An advantage of basing the verification of the theoretical results with experiment on the 
j7'(2)_]2ioments is that this quantity does not depend on the exit spin (the spin corresponding 
to the lowest observed transition in an SD band) while many other quantities do. For this 
reason we intend to focus on the rotational frequency rather than spin as an independent 
variable in function of which several microscopically calculated quantities will be displayed. 
However, the question of the relative spins in the multi-band problem, or equivalently, of 
the relations between the exit spins in the bands observed, is one of the most urgent at 
present. For this reason a special section will be devoted to discussing the total energy vs. 
spin dependence, £ = £{I). 

Let us mention here that, similarly to numerous other studies devoted to the inter- 
pretation of the intrinsic configurations of the SD bands, our work will not be able to give 
definite answers in several points. We feel, nevertheless, that our method allows to pin down 
relatively reliably at least the possible alternative interpretations of the bands in question. 

While the lack of knowledge of the exit spins remains one of the most serious experimental 
problem, many analyses of the SD nuclei have benefited from the strong similarities between 



various bands in various SD nuclei (the "identical band" correlations, cf. Ref. ||T3| for a 
recent review) . By comparing the theoretically predicted structures of the bands related by 
the "identicality" (or "sameness") relations, it was possible to confirm or eliminate certain 
hypotheses related to the underlying single-nucleonic structure. This type of correlations 
are abundant in the ^^^Tb nucleus and will be used to further test the microscopic results 
and predictions. 

It is expected that some amount of pairing correlations remains in the SD bands up to the 
highest spins measured and that the pairing correlations decrease with increasing angular 
momentum. Although, as we will demonstrate it, in the ^^^Tb nucleus the most important 
properties of the SD bands such as, e.g., the behavior of the dynamical moments in function 
of the rotational frequency can be described reasonably well without taking into account 
the pairing, the systematic small discrepancies which exist, are increasing when the pairing 
correlations are taken into account. Moreover, the arguments will be given demonstrating 
that at least in the case of some parametrizations of the nuclear interactions another source 
of small systematic deviations between theory and experiment may be attributed to slightly 
exaggerated nuclear volumes as calculated within the mean field approaches (cf. Sec.|V| for 
details). 

A number of structural features related to the high-spin behavior of the SD nuclei in 
the A~150 region have been discussed in the literature. Early studies were based on the 
Nilsson or WS cranking approaches, and performed without pairing WM or with the dynamic 
pairing correlation [p!5|-[T8[l taken into account by using either the random phase approxima- 
tion (RPA) or the particle number projection (PNP) together with the monopole pairing 
interaction. In most of those early studies the experimental information has been strongly 
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limited, as compared to the amount of data available today. In more recent studies, the 
WS |]l9l or the Hartree-Fock-Bogolyubov (HFB) methods have been applied to a few 
low-lying bands in a number of selected A~150 nuclei with pairing correlations treated in 
terms of the Lipkin-Nogami method. 

The present study differs from the previously published ones in the following: (i) due 
to a rich experimental information available now, we can concentrate on systematic dis- 
crepancies between the theoretical and experimental data, (ii) we discuss the analogies and 
differences between the self-consistent (HF) and the non-self-consistent (WS) cranking ap- 
proaches within the same analysis (applied simultaneously on a common set of experimental 
data), (iii) we analyze simultaneously all available data in the nucleus in question, (iv) by 
applying the particle-number-projected pairing approach we are able to show that the sys- 
tematic discrepancies cannot be attributed to the pairing correlations, and consequently (v) 
we arrive at a suggestion that a modification of the underlying interactions will be needed 
in order to improve the agreement with data. 



II. SUMMARY OF THE EXPERIMENTAL RESULTS 

As mentioned earlier, in the ^^^Tb nucleus eight SD bands have been observed. In ad- 
dition to the two SD bands discovered relatively early pO| , pT| , six new SD bands have been 



found and assigned to this nucleus recently p2| , |23[| using the EUROGAM spectrometer. 



Those recent data were obtained from experiments at the Nuclear Structure Facility, Dares- 
bury Laboratory, using the ^^°Te(^^Al,6n)^^^Tb reaction at a bombarding energy of 154 MeV. 
The details of the experimental procedure can be found in P^ , |2B[] ; here we only recall that 



the 7-rays were detected by the EUROGAM spectrometer which consisted, for this 

experiment, of 42 operational Ge detectors. The data rates in the individual Ge detectors 
were 7-9 kHz. The trigger condition used to select unsuppressed events was that at least 
7 Ge detectors were involved in a coincidence. After the suppression requirement was in- 
cluded, the average number of coincident detectors dropped to ~3.9. A total data set of 
5.5x10^ suppressed events with Ge fold > 3 was recorded and after unpacking the high- 
fold-multiplicity events led to 3.1 xlO^ 7^ coincidences. The final spectra were obtained by 
summing combinations of 3-dimensional gates set on the SD transitions in a 4-dimensional 
analysis ||2^ . 



The assignment of the new 6 excited SD bands to ^^^Tb was based on an unambiguous ev- 
idence of known 7-ray transitions of the normal deformed nucleus following the de-excitation 
from the SD band members. The SD bands have been labeled 1,2,3,. . .,8. The previously 
reported SD bands have been labeled 1 (called the yrast band, although its position in en- 
ergy with respect to other bands is not experimentally established) and 2 (the band which 
shows a remarkable similarity in terms of the 7-ray energies with the SD Band 1 of ^^^Dy). 
The other excited SD bands have been arbitrarily attributed labels ranging from 3 to 8. 

Employing the same reaction, the ^^^Tb nucleus has been recently reinvestigated us- 
ing the Eurogam II array installed at the Vivitron accelerator of the Centre de Recherches 
Nucleaires in Strasbourg. This array consisted of 54 high-efficiency Ge detectors (30 ta- 
pered coaxial detectors and 24 clover detectors) providing a photopeak efficiency of 7.4% at 
1.33 MeV. With the large statistics and high peak-to-background ratio, it has been possible 
to obtain a higher precision of the SD 7-ray transition energies which in the present study are 
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used to extract accurate values of the corresponding dynamic moments. The experimental 
results for Bands 1 and 2 are presented in Table |. 

The intensities (with uncertainties) of the different bands relative to the yrast band are: 
50(5)% (Band 2), 35(5)% (Band 3), 6(2)% (Band 4), 10(2)% (Band 5), 9(2)% (Band 6), 
11(2)% (Band 7), 7(3)% (Band 8). 

Let us also mention that in the determination of the experimental frequencies and dy- 
namic moments we use the expressions 
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_ E^{I + 2) + E^{I) 



(2.1a) 

^ ' E{I + 2) -2E{I) + E{I -2) 

(2.1b) 



E^{I + 2) - E.,{I) 
where E^{I)=E{I)-E{I - 2). 



III. WOODS-SAXON CALCULATIONS 



In this Section we summarize the main points of the calculations based on the average 
nuclear field with the deformed WS potential. The following three subsections present: 
(i) the structure of the Hamiltonian used, Sec. |III A| , (ii) some remarks related to the self- 
consistency in the case of a priori non-self-consistent formalism with the WS potential, 
Sec. [1IIB| , and (iii) a short description of the actual variant of the Strutinsky method applied. 
Sec. [Ill (J| The details can be found in the literature cited below and will not be repeated 
here. Even though most of the mathematical expressions are available in the literature, there 
exist differences in notation between various authors, and we found important to recall the 
definitions of some basic quantities. 



A. Deformed Woods-Saxon Hamiltonian and the cranking approximation 



The single- nucleonic energies and wave functions for nucleons moving in the deformed 
potential of the WS type, 

V{r; Vo, n, a, ro) = ^ ° ^ (3.1) 
{1 + exp[distso(T';ro)/a]} 

are found by solving the corresponding Schrodinger equation. The spin-orbit term is included 
in the usual deformation-dependent form: 

2 



h 



2mc 



X 



Vo[l±K 



N-Z l 
N+Zi 



{1 + exp[distsso(^^; '"so)/aso]} 



xp} -s, 



(3.2) 
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and the Coulomb potential is added for protons. In the above relations the "+" sign applies 
for protons and the "— " sign for neutrons, while p and s denote the linear momentum and 
the spin operators, respectively. For reasons of compatibility with previous studies, we keep 
the factor [h/{2mc)]^ in the definition of the spin-orbit potential (m is the nucleonic mass). 
Then the spin-orbit strength parameter. A, is a dimensionless quantity. We use A = 35 for 
neutrons and A = 36 for protons. 

The nuclear surfaces Sq and Sso, defining the central and spin-orbit fields, respectively, 
are expressed in terms of the spherical harmonics as 



^0{SO) 



R{9, 



0(30) 



c a 



1 + H "Am^Am( 
A=0 



(3.3) 



In the above relation, c{a) is a function of all the deformation parameters that guaran- 
tees the constant volume for any deformation, and R, 



0(SO) 



ro(so)^"'^^^ with the central- 
potential radius parameter rg = 1.347 fm for neutrons and Tq = 1.275 fm for protons; sim- 
ilarly rso = 1.310 fm for neutrons and rso = 1.200 fm for protons. The central-potential 
depth-parameters are: Vq = —49.6 MeV and k = 0.86 while the diffuseness parameters are 
a = 0.70 fm for the central part and a^o = 0.70 fm for the spin-obit part of the potential, 
independently of the isospin (for details see Refs. and [0]). Finally, the geometrical 



distance of a point r from the auxiliary surfaces So(so) defined above is denoted by dists^^gQ^ ; 
the latter quantity is negative for r in the nuclear interior and positive otherwise. 

In our calculations we allow for the total Routhian (for definition of the Routhians see 
the following sections) minimization with the deformation-mesh defined by 



and 



aao = 0.400, 0.425, . . . , 0.800 (17 points) 



040 = 0.020, 0.040, . . . , 0.180 ( 9 points) 



(3.4) 



(3.5) 



with corresponding to 5 "planes" for fixed 7 deformations [0:20 = /5cos(7) and 0:22 = 
P sin(7)] as follows 



7 = -10°, 5°, ...,10° (5 points). 
We solve the cranking equations of the standard form 

for 



hu = 0.00, 0.05, . . . , 1.50 MeV (31 points) 



with 



T + V + Vso + Vc - ujI, 



(3.6) 
(3.7) 
(3.8) 
(3.9) 



where the WS potential V and the related spin-orbit term Vso are defined in Eqs. (^TTj) 
and ( |3.2|) while the Coulomb potential Vq is that of the classical uniform sharp-edge charge 
distribution corresponding to the actual nuclear shape. 
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As always in numerical calculations of the matrix elements used for constructing the 
Hamiltonian matrix, both the number of basis states and the number of points in the Gauss- 
Hermite integration formulae have been selected in such a way that the final results do not 
depend in any significant way of that choice (for more details see Ref. 0). In our case the 
number of the deformed harmonic-oscillator states used in each of the two parities was about 
300. 



B. Self-consistency between the potential and density 



In the following, the results of the cranking calculations performed within the HF ap- 
proach will be compared with those obtained for the deformed WS potential. Although 
it is not possible to formulate the self-consistency concept in the two approaches in the 
same way, and the WS calculations are a priori considered non-self-consistent, it will be of 
some importance to recall that our way of defining the deformed WS potential imposes a 
correlation between the equipotential surfaces 



V{r) = const., 
and the surfaces of the constant density defined by 

Pi^) = J2^nir)^n{r) = const. 



(3.10) 



(3.11) 



as discussed in detail in . The second-order Thomas- Fermi expression |^9|] for the density 
functional reads 



p(r) 



1 



37r2 



2m 



;2 [-f 



3/2 



(3.12) 



8 2m [ 4[X-V{r)f 



+ 



AV{r) 



[X - V{r)y 



where the terms of the order 0{h'^) have been neglected. Requiring that the potential is 
equal to a constant, cf. Eq. ( |3.10D , implies by virtue of Eq. ( |3.1| ) that 



Using the fact that 



and that 



dist(r) = const.. 



Vdist(r) • Vdist(r) 



(3.13) 



(3.14) 



Adist(r) = const. (3.15) 

for r belonging to the equipotential surface, one obtains after elementary transformations, 
cf. Ref. [^], that also 

W{r) ■ VV{r) = const, and AV{r) = const. (3.16) 
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As a result, at least within the accuracy of 0{h'^) in the Thomas- Fermi approximation, 
condition ( ^.1(J| ) implies condition ( |3.11| ). Such a correlation between the geometry of the 
potential and that of the density is not valid generally. Neither the Nilsson model nor a 
deformed WS potential with another type of argument in the exponent of Eq. (^.11), as used 
sometimes by other authors, gives such a correlation. 

Obviously, the spin-orbit term may modify the above result (cf. e.g. Ref. [0), but since 
this term has, on the average, matrix elements which are nearly an order of magnitude 
weaker as compared to those of the central part, its effect in the discussed context is also 
expected to be relatively weak. 

Anticipating some of our results presented below, and in particular surprising similarities 
between the high-spin behavior properties obtained by using the self-consistent HF method 
and a non-self-consistent approach with the deformed WS Hamiltonian, one should bear in 
mind the particular property of definition ( |3.1D used in this study. 



C. Strutinsky method and the total-Routhian calculations 

There are numerous presentations of the Strutinsky method, as given for example in 
the original articles and then applied to finite potentials (as e.g. in Ref. Here 



we closely follow the presentation given in Ref. [g], and we only briefly discuss the total- 
Routhian-surface calculations, for which several authors use prescriptions slightly differing 
in details. 

Our calculations of the total nuclear Routhians are based on the expression 

pa; _ rpuj=0 , rpu)=0 , \ ^ f u) ui=0\ fn i 

^ -'^macro -'^micro / , V^n )i \o.L I ) 

n 

where the summation extends over all the occupied nucleonic states (i.e., in particular takes 
into account the particle-hole excitations with respect to the reference configuration intro- 
duced in Sec.[V|). The microscopic (or shell) energy is defined as usual by 

E-^L = E - < E >.heu (3-18) 

where the summation over u extends now over all the single-nucleonic states of the reference 
configuration. The last term denotes the results of the Strutinsky averaging procedure. 
This averaging depends on two standard parameters: the smoothing-polynomial order p 
and the averaging width 7av. In our applications we have p = 6 and 7av = 3.2 husheii with 
^i^sheii = 41/^4^/^ MeV. Arguments related to the stability of the final result with respect to 
(small) uncertainties in these parameters can be found in several publications, see e.g. Ref. 
101 and references therein. 



IV. HARTREE-FOCK CALCULATIONS 

In the following Sections, some details of a particular realization of the HF approach in 
our study will be described; most of the formal development can be found in the literature 
and are also summarized in Ref. |Q . Section [IV A| presents some technical aspects related 
to the harmonic oscillator basis employed, the numerical calculations of the jT^^^-moments 
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are discussed in Sec. |IV B| , and the use of the quadrupole constraint to stabihze the iteration 
process is briefly commented in Sec. IV C . 



A. Harmonic-oscillator basis 



In the present study we have performed the HF calculations of SD rotational bands 
using the Skyrme effective interaction. For this purpose a new numerical code HFODD 



has been constructed |3J] which uses the three-dimensional Cartesian deformed harmonic 
oscillator (HO) basis. The HO basis has up to now been used in solving the HF equations 
with the Gogny interaction p5|-[39[]. On the other hand, for the Skyrme interaction this 



method has been used at early stages of investigations and then abandoned in favor of 
the calculations in spatial coordinates 

The main advantage of using the spatial coordinates consists in the fact that one may 
use the same spatial grid of points to treat nuclei of different deformations, while in the 
former method the basis parameters may need to be adjusted with varying deformation. 
(The sensitivity of the results with respect to this kind of a dependence decreases with 
an increasing harmonic-oscillator basis cut-off; in principle the problem disappears in the 
limit of an infinite basis). However, the SD states have very similar deformations in many 
neighboring nuclei, and in such a case this particular deficiency disappears. We may then 
profit from the fact that the calculations which use the HO basis are in general much more 
stable, have better convergence properties, and hence require less computational effort. 

The details concerning the HFODD code will be presented elsewhere 0]; here we only 
give a few of its basic characteristics pertaining to the present application. The calcula- 
tions have been performed using a fixed basis given by the HO frequencies huj ±=11. 200 MeV 
and ^a;|| =6.246 MeV in the directions perpendicular and parallel to the harmonic-oscillator 
symmetry- axis, respectively. These values have been obtained by standard prescriptions 
developed for diagonalizing the deformed WS Hamiltonian in the HO basis, and corre- 
spond to the WS potential with the deformations /?2=0.61 and ^4=0.10. The basis has been 
restricted to the fixed number M of HO states, i.e., keeping the lowest M single-particle 
energies eHo={nx+''T'y+l)huj± + {nz+^)huj\\. The actual calculations have been performed 
with M=306. This corresponds to the maximal numbers of oscillator quanta equal to 8 and 
15 in the perpendicular and parallel directions, respectively. 

The stability of results with respect to increasing the size of the HO basis is the main 
concern in all analyses using this technique, and have been studied in detail in Ref. 



In fact it is well known that the total energies are very slowly convergent as functions of 
the basis size. Indeed, in the present calculations a similar effect is seen when we increase 
the basis from M=306 to M=604, which corresponds to reaching the maximal numbers 
of 11 and 20 HO quanta in the two main directions. Then the total energy of the yrast 
configuration of ^^^Tb at the rotational frequency of huj=O.Q MeV decreases from —1199.322 
to —1201.889 MeV. A similar mechanism is well known also from other type of studies, with 
the Schrodinger equation diagonalized in a given basis (as for instance with the WS or the 
folded Yukawa- type potentials). Fortunately the relative (excitation) energies become stable 
for smaller bases than the absolute values . 



None of the interesting rotational properties depends in a significant way on an increase 
of the basis from M=306 to M=604. For example, we obtain the changes of the average 
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angular momentum from 59.95 to 59.99^, of the mass quadrupole moment from 39.82 to 
39.90 b, and of the mass hexadecapole moment [41 1 from 4.21 to 4.24 b^. The relative changes 
as function of hu are even smaller. Most interestingly, the JT"*-^-* -moment changes at the same 
time only from 84.71 to 84.84 MeV'^ In view of the fact that the calculations for M=604 
are a factor of five slower than those for M=306, we consider it reasonable to keep the latter 
size of the basis in all calculations presented. 



B. Calculations of the dynamical moments 

An important advantage of using the HO basis consists in the fact that a rather small 
number of iterations is sufficient to obtain a high level of convergence. For example, in the 
case of the yrast configuration in ^^^Tb at huj=0.6MeV, starting from the rotating WS state 
and performing 50 iterations is enough to obtain the total energy accurate up to 0.5 keV. 
For all values of rotational frequency and with the same number of iterations performed, 
the relative precision remains comparable to the one cited above. 

This is very important for the calculation of the second moment of inertia J'^'^\ This 
quantity has to be calculated by a numerical differentiation, either of the average angular 
momentum I W3, 



/ = (^"^141$^), (4.1) 
with respect to the rotational frequency u, Eq. ( |1.1|) , or of the total energy S, 

£ = (4.2) 



with respect to u: 



^ =-7r + (^^2o) - Q20) (4.3) 
uj auj ^ ^ uj auj 



where the cranking wave functions |\1''^) are obtained from the HF variational equation 



5{^-\H + C2 ((Q20) - 020)' - = 0- (4.4) 

In order to obtain accurate values for J''^'^\ the relative errors of I and 8 as functions of uj 
should be kept as small as possible. The difference of values of J^"^^ calculated from Eqs. ( |1.1D 
and ( [4.3| ) is a measure of self-consistency of solutions. In our example, the previously quoted 
value of iT'*^^'' =84.71 /i^ MeV~^ has been obtained from Eq. ( [1.1|) , while at the same time 
Eq. ( [4.3|) gives J^(2)=84.66^^MeV-^ These resuhs have been obtained from the first-order 
finite-difference approximation of derivatives between ^a;=0.60 and 0.61 MeV. Together with 
the discussed uncertainties related to the restricted size of the basis, these values indicate that 
our calculated dynamical moments are accurate to about O.l^^MeV"^, which is generally 
a higher precision than that presently available for the experimental values of JT'-^-* . 
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C. Quadrupole constraint 



The SD states are localized in the pocket of the potential energy surface which is several 
MeV above the main minimum. When the HF equations are solved by the standard itera- 



tion procedure which we also use in our calculations, the wave functions may become 
gradually polluted by components of the solutions residing principally in the main well, and 
this may slow-down or even preclude any convergence to the secondary minimum. In order 
to prevent this kind of effects we have included in the minimized functional a quadratic 
constraint on the quadrupole moment [Q, i.e., we minimize 



£' = ^ + c2((g2o)-g2o)' (4.5) 

with a small stiffness parameter 6*2=0.01 MeV/b^ and the constraint value of Q2o=42b. 
This is a very weak constraint which has no visible effect on the self-consistent values of the 
quadrupole moment. It shifts, however, the normal-deformed-minimum energy sufficiently 
high, so that it does not present any undesired competition with the SD one. 

The quadrupole constraint is at the origin of the second term in Eq. ( [4.3| ) , which compen- 
sates for the fact that the quadratic constraint adds to the functional a term which depends 
on the rotational frequency. The magnitude of this second term depends on the variations 
of the average quadrupole moment 

(Q20) = (^1Q2o|^") (4.6) 

along the rotational band (see Sec.0) and on the difference between {Q20) and the tar- 
get value (520- In the present calculations this term can be for some bands as large as 
0.5 /i^ MeV~^, and, therefore, cannot be neglected when comparing the values of JT'-^-* calcu- 
lated from Eqs. (g^ and (p|). 



V. THEORETICAL RESULTS AND COMPARISON WITH EXPERIMENT 

In this section we present the results of the HF calculations of the multi-band excitation 
structure of the ^^^Tb nucleus. In several cases the cross-check results based on the deformed 
potential of the WS form will also be presented. We divide this discussion into several steps: 
Sees. |V A| and V Ij are devoted to the analysis of the rotational bands in terms of the single- 



particle Routhians, total energies, and dynamical moments J''^'^\ while the next ones present 
problems related to exit spins, shape changes along the rotational bands, and pairing. 



A. Superdeformed yrast band in Tb 



Most probably, only the lowest energy SD bands are populated in the present-day exper- 
iments on the SD nuclei. More precisely, the bands lying not too high in the energy scale, 
at spins between, roughly, 60 Ti and 80 Ti, can actually be detected. Unfortunately not much 
more can be said about the population selectivity and no rules for, for instance, the pop- 
ulation probability as a function of the actual configuration have been established. Under 
these circumstances the only way of performing a systematic theoretical analysis consists 
in calculating numerous bands in excess of the number of those actually measured and in 
proceeding with a more detailed and critical selection analysis next. 
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1. Single-nucleonic energy spectra 



As the first step the single-particle spectra of the deformed WS cranking Hamiltonian, 



cf. Sec. |TJ, have been calculated as functions of the cranking frequency for hu ranging from 
to IMeV for a few points corresponding to the typical quadrupole 020 ~0.58 and the 
hexadecapole 040 ~0.10 deformations known to represent well, on the average, the SD 
nuclei in the y4~150 mass range. Knowing those levels, we define the reference configuration 
that corresponds to the lowest-level occupation scheme for possibly largest fiuj range. The 
symbols of the particle-hole configurations used below refer to such a specific reference 
configuration as the underlying structure. It is worth mentioning that often more than one 
configuration could quite well serve as the reference and that the actual choice of such a 
configuration is a matter of commodity. In fact, an inconvenient choice of the reference 
would merely result in a more complicated labeling of the excited bands: a 1-particle 1-hole 
excitation with respect to a given reference configuration may become 2-particles 2-holes, 
3-particles 3-holes etc. excitation with respect to the inconveniently chosen reference. 

In the present analysis we consider Hamiltonians having two discrete symmetries: the 
parity ^ and the y-signature Ry 

Ry = exp^—iniy), (5.1) 

which has the eigenvalues r=+i or r=—i. The reference of our choice is defined in terms of 

four types of the occupation numbers: A^++, N^_, iV__|_ and referring to positive-parity 

and signature r=+i, positive-parity and signature r=—i, etc., respectively. We have two 
sets of such numbers for protons and neutrons each. In ^^^Tb the self-consistent solution for 
the reference configuration is defined by occupying the single nucleonic orbitals according 
to the following scheme 



(5.2) 



In Fig. we show the single-particle Routhians for protons and neutrons obtained with the 
SkM* parametrization of the Skyrme interaction for the reference configuration given in 
Eq. ( [5.21 ). All HF calculations are in the present study performed using the complete Skyrme 
functional, cf. discussion in Ref. 0. 

Let us remark that the interpretation of the experimental SD ^^^Tb(l) yrast band in 
terms of the configuration implied by relations ( |5.2| ) corresponds to the occupation of all 
the lowest-lying proton Routhians with the last orbital occupied 7r[6,5,l]3/2(r=— i) (at large 
uj taking another label: 7r[6,4,2]5/2(r=— i), cf. Fig.|I|). This reference configuration has the 
lowest calculated total energy at a fixed spin, and therefore, from the theoretical point 
of view, we can call it the yrast SD configuration in this nucleus. Configuration assign- 
ment (|5.2|) is confirmed by experimentally observed band similarities between ^^^Tb(l) and 
^^°Gd(2), whose several transitions correspond as being nearly identical. Indeed, this obser- 
vation can be viewed as a direct consequence of exciting the last occupied proton orbital, 
7r[3,0,l]l/2(r=+i) in ^^'^Gd{2), into 7r[6,5,l]3/2(r=— i), the last occupied proton level in 



Ar++ 


= 22, 


N+- 


= 22, 


N-+ 


= 21, 


A^„_ 


= 21, 


N++ 


= 15, 


iV+- 


= 16, 




= 17, 


A^__ 


= 17. 
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^^^Tb(l). Since the alignment properties and the jT'^^^-contributions originating from the 
7r[3,0,l]l/2(r=+z) orbital are known to be negligible, Refs. , the assignments of 

both configurations in the two nuclei can be viewed as confirmed. 

In order to establish the uncertainty margin due to various parametrizations of the 
interactions that have been used in the literature, the results analogous to those in Fig.|l| 
but for the SIII parametrization of Ref. ^8|] are presented in Fig. ^ One can see that the 
characteristic gap structures for the protons (in the form of the sequence of gaps at Z = 66, 
Z = 64, and Z = Q2 visible in Figs.|l| and ^ are very similar in both parametrizations. 
These gap structures are defined by very similar positions of the 7r[6, 5, l]3/2 and 7r[3, 0, l]l/2 
orbitals. Other proton orbitals, lying above and below those gaps, differ only slightly in the 
level ordering. The differences at u; = in the level positions do not exceed typically ~ 500 
keV. 

It will be instructive at this point to compare the single-particle spectra obtained by 
using the HF method with the analogous results obtained by making use of the deformed 
WS cranking Hamiltonian. For this purpose the three-dimensional total-Routhian surfaces 
for the ^^^Tb have been calculated in function of the quadrupole (/3,7) and hexadecapole {P^) 
deformations, cf. Sec. |III A| . The total-Routhian surfaces have been calculated in function of 
the rotational frequency and the minimum-deformation path composed of sets of (/3,7,/?4)- 
deformation parameters that vary with increasing huj has been found. The deformed WS 
cranking single-particle spectra obtained in this way are presented in Fig.|^. 

Let us observe by comparing Figs.|l|, ^, and |^, that the global features of the single- 
nucleon spectra obtained by making use of the deformed average field model with the WS 
Hamiltonian and those of the HF with the SkM* and SIII interactions are quite similar. 
What is perhaps most striking is the fact that the WS Hamiltonian and the HF method with 
the SIII parametrization give very similar single-nucleonic levels indeed. The similarities are 
particularly well visible in the neutron spectra in Figs. |I] and ^, and only to a slightly lesser 
extent in the proton spectra displayed in the Figures. This observation can be interpreted 
in favor of both techniques that employ the two different average fields thus increasing our 
confidence in the implied results of the microscopic calculations. 



-moments 



Since the exit spins of most SD bands are not known, at present the only possible way of 
testing the calculated S versus / relation by experiment consists in comparing the resulting 
j7'(2)_]2ioments. Before proceeding systematically with such comparisons some additional 
comments will be useful. We begin by discussing the WS calculations which have been 
performed here using the so-called "universal" parameters. The corresponding parameter 
set has been introduced in Ref. [^, and later tested in many detailed calculations of the 
high-spin properties in nuclei. The corresponding parametrization of the central part of the 
WS potential employs the radius-parameter values of ro(z/) = 1.347 fm and ro(7r) = 1.28 fm 
for the neutrons and protons, respectively. Historically, these parameter values have been 
fitted to improve the quality of the single-particle level spectra rather than to reproduce 
well the geometrical properties of nuclei, and it happens that they exceed considerably the 
more commonly accepted "geometrical" value of ro — 1.24fm. 
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One of the direct implications of the above difference is that the kinematical moments 
of inertia 



J« = -, (5.3) 

UJ 

calculated from, e.g., the Fermi-gas model with the geometry defined by using the above 
radius-parameter values systematically exceed the rigid-body moments of inertia calculated 
from the classical formula with the standard nuclear matter density. Within a classical but 
also within a semiclassical model, the moments of inertia simply scale with the square of the 
radius parameter, Tq. As a result one should expect that the moments of inertia in various 
bands, calculated by using the "universal" parametrization of the WS Hamiltonian will be, 
on the average, too large. We can see that such an expectation, valid for the kinematical 
moments will apply directly to the dynamical moments ( |1 . 1|) , 

J^'^ = ^ = J'~'^+J-^, (5.4) 
duo auj 

and a scaling of the J' -moments implies a similar scaling of the J'^'^^ ones. 

By scaling the calculated J'^^^ and JT"*-^^ -moments we should, as a consequence, also scale 
the related physical quantities: rotational frequency, energies and Routhians. One can see 
this after noticing that the scaling of the JT"*-^^ moment with a factor /, that satisfies relation 
( p.3[ ), for the quantized angular momenta J = 0, 2, 4, . . . (which remain independent of any 
scaling), implies the necessity of scaling the rotational (cranking) frequency by a factor of 
/^^. It then follows that the energies and Routhians also scale as f~^, so as to consistently 
satisfy the canonical relation E = R + hul. 

The results of the calculations of the JT*^^^ moment for the reference (in this case coinciding 
with the yrast) configuration of ^^^Tb are given in Fig. |, where the WS and HF curves are 
compared with experimental data. Top part of the Figure gives results with no scaling 
applied. One can see, that the calculated values systematically exceed the experimental 
ones. In the WS case we know that, as mentioned above, the numerical values of the radius- 
parameters are too large and we can clearly consider at least a part of the discrepancy as 
related to "too large a volume of the nucleus". Therefore, in the bottom part of Fig. ^ we 
show the same jT^^^-moments scaled by the factor /=0.9. 

In the HF case there is no obvious discrepancy between the calculated and experimental 
nuclear radii. Nevertheless, the calculated J' ^"^^ -moments as well as quadrupole moments are 
systematically slightly too large (see Refs. [^,0 for recent measurements of the quadrupole 
moments and Ref. for the results of systematic HF calculations of the quadrupole mo- 
ments and comparison with experiment). Such an overestimation is consistently obtained for 
many different Skyrme force parameters f^. At present, no explanation of this systematic 
effect can be found. For the sake of the presentation we have therefore decided to use the 
same phenomenological scaling factor /=0.9 to present the results of HF calculations. 

This scaling procedure, as we have checked by an explicit comparison with the experi- 
mental results, works well for all the bands in the ^^^Tb nucleus simultaneously. Moreover, 
we have checked that by using the same scaling factor /=0.9 for many bands and for several 
nuclei simultaneously one reproduces well the experimental JT*^^-* moments which otherwise 
remain too large. This conclusion is based on comparisons for the yrast bands in ^^^Dy, 
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^^^Dy, ^^^Dy, ^ Tb, and ^^'^Tb, and also for a few dozens of excited SD bands in the mass 
y4~150 region. Light Gadohnum nuclei (^^^Gd and ^'^''Gd) have strong alignment effects in 
their yrast SD bands and are not well suited for the tests discussed here. Slopes of the J^'^^ 
moments in ^^^Gd, ^^^Gd, and ^^''Gd are overestimated in the calculations; for the discussion 
of the latter nucleus see Sec |V E| . Conversely, not introducing the aforementioned scaling we 
find the calculated values systematically too high, by 8 to 10%. 

After applying the scaling procedure according to the discussion above, (cf. Fig.^ for 
illustration), the jT^^^^-moments are multiplied by the factor of / and the angular frequen- 
cies by a factor of /~^. Then, a very good correspondence with the lo dependence of the 
experimental data is worth noticing. Below we present all results scaled according to the 
prescription and parameters defined above. Let us mention that no physical conclusion of 
this paper depends on using or not using the scaling procedure described above; the scaling 
can be viewed as merely a "graphical" trick to facilitate a comparison with the experimental 
data. Moreover, we have not attempted any detailed fit of the scaling factor /=0.9 which 
can still be varied within about ±0.02 without any obvious general deterioration of the 
agreement with data. 



B. Excited SD bands in ^^^Tb 



The problem of the definition of the yrast configuration discussed in Sec. |VA 1| did not 
pose any serious difficulty. Firstly, in nuclei with the nucleon numbers in the vicinity of 
large SD gaps (as e.g. Z=66 and A^=86) the yrast configurations corresponding to the 
occupation of the lowest-lying levels are relatively easily and unambiguously determined 
over a relatively long cu-range since no, or very few level crossings occur at the Fermi level 
(cf. Fig. |1]). Secondly, the experimentally dominating SD band (i.e., the one with the highest 
population intensity) is also easily distinguished. In the case of the other SD bands neither 
of the observations applies. In experiment several bands are populated with comparable 
(usually weak) intensities and their relative excitation energies are totally unknown. In 
theory, numerous close lying bands remain a 'priori the candidates for a much smaller number 
of the experimental bands to be associated with. 

Moreover, the uncertainties in the relative positions of the theoretically calculated levels 
which sometimes influence the level ordering may reach 500 keV or more. This estimate 
corresponds to a typical discrepancy between the results obtained within commonly used 
average field models (like WS potential or Skyrme interaction parametrizations, SkM* and 
SIII as discussed above and illustrated in Figs. |], || and |^) and may, occasionally, become 
markedly larger. Consequently, in analyzing the structure of the nucleus in question we have 
applied the following approach: 

(i) Our primary goal in this section is to recognize within the applied formalism the 
most probable single-nucleonic (particle-hole) structure that underlies the behavior of the 
experimentally found bands. 

(ii) We are using, for the sake of qualitative argumentation, the calculated ordering of 
various SD bands: the lower the band energy at the high-spin limit, the higher the chance 
to find such a band among those populated in experiment. 

(iii) We use the experimental j7^*^^^-moments and verify the possible single-nucleonic con- 
figurations by directly comparing J'^'^\uj) with experiment. 
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Another important piece of information suitable for identifying single-particle configu- 
rations is the sameness of bands in different nuclei, i.e., the equality of the corresponding 
dynamical moments and specific relations which occur for the relative alignments, see the 



review in Ref . ||T3[ . In the present work we mainly use and analyze the properties of the dy- 
namical moments, because those of the relative alignments are very sensitive to the details 
of the effective interactions p|. A readjustment of the interaction parameters, especially 
those related to time-odd components of the mean field, seems to be necessary for a precise 
theoretical determination of the relative alignments. 

In order to apply the guidelines (ii) and (iii) above, we first present the results for the 
lowest-lying one-particle one-hole excited bands. In Figs. ^ and ^ we show the total energies 
for the proton and neutron excitations, respectively. Here and in the following we label the 
configurations by the corresponding Nilsson labels of hole and particle states, as well as by 
the total parity and signature quantum numbers supplemented by (arbitrary) consecutive 
number attached to each configuration. Since in ^^^Tb the size of the proton shell gap is 
smaller than that for the neutrons, there are more low-lying proton than neutron excited 
bands. Since we know that there exist uncertainties in the single-particle level positions 
that may influence the relative positions of some bands in a non-negligible manner, the 
results in Figs. ^ and |^ should be treated as a semi-quantitative guide. The Figures show 10 
proton particle-hole and 12 neutron particle-hole excited bands, respectively, that appear 
to be the lowest at / ~ 80^, i.e., in the spin range where the superdeformed bands are 
populated. We have placed some of the curves in the upper parts, and some others in 
the lower parts of the Figures trying to illustrate an existence of "families" of bands, each 
family characterized by some specific features like slopes or band crossings. In the case of 
neutrons, the separation into two families is particularly well visible. The differences in this 
case (cf. also Figs. ^ and |^) correspond to the structure differences in the two characteristic 
hole-orbitals (z/[6, 5, l]3/2(r=+z) or z/[6, 5, l]3/2(r=— i), top parts, and z/[7, 6, l]3/2(r=— i) 
bottom parts). 

Figures |^ and ^ present the theoretical predictions for the dynamical moments corre- 
sponding to the bands shown in Figs. |^ and ^, respectively. One can see that the theoret- 
ically calculated moments form "families" characterized by a similar behavior also in the 
proton case. In particular the bands built on the configurations of the type 7r[3, 0, l](r=±z) 

7r[6, 6, 0](r=+i) and 7r[6, 5, l]3/2(r=— i) — > 7r[6, 6, 0](r=+z) are characterized by similar 
slopes, in terms of the J'^'^^ moments; these slopes are different from that of the refer- 
ence band. On the neutron side, the similarities between the bands related to exciting 
the iy[7, 6, l]3/2(r=-i) orbital into the orbitals i/[5, 1, 4]9/2(r=±i), z/[5, 2, l]3/2(r=±i), and 
i/[4, 0, 2]5/2(r=±i) deserve noticing. Distinguishing among these properties will be partic- 
ularly useful when attributing the configurations to the experimental bands since, as men- 
tioned before, several theoretical configurations provide similar J'^'^^ moment dependence on 
the rotational frequency, while their energies are also quite close. 



1. Bands ^^^^^{2), ^^^Tb{3), and ^^'^Tb{4) 

We proceed with the interpretation of the experimental results for Bands 2, 3, and 4. 
Although the most likely structures of these bands have already been discussed before in the 
literature [iB,23|,P,r^|, we found it instructive to illustrate the procedure of attributing the 
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theoretical bands to the experimental ones, the procedure which will be of more importance 
for other bands for which an adequate interpretation remains less easy. 

By eliminating the bands whose i/'^^^-moments differ relatively strongly with respect 
to the experimental result for Bands 2-4, we arrive at a comparison presented in Fig.|^. 
Here and in the following we denote the HF states by their dominant Nilsson labels at 
high frequencies, cf. Fig.|l]. For example, the above mentioned 7r[6, 6, 0]l/2(r=+z) Routhian 
originates at low frequencies from the vr[6,5, l]3/2(r=+i) state. 

Experimental results for Bands 2 and 4 agree rather well with calculations for config- 
urations based on exciting either of the signature partners of the 7r[3, 0, l]l/2(r=±i) or- 
bitals into the 7r[6, 6, 0]l/2(r=+i) state. As seen in Fig. ^ the configuration based on the 
7r[3, 0, l]l/2(r=+i) Routhian has the energy lower by about 1 MeV than that based on the 
corresponding signature partner, and therefore the former should be retained as the most 
likely interpretation of Band 2. Similarly, the dynamical moment of the 7r[6, 5, l]3/2(r=— i) 

7r[6, 6, 0]l/2(r=-|-i) configuration, Fig.|^, middle part, reproduces very well the data for 
Band 3, and this supports the corresponding interpretation (proposed also in Ref. [|23|). The 



suggestion of this configuration is based simultaneously on the fact that it corresponds to a 
low-energy excitation, as it is seen from Fig. |^, and that it provides a good agreement of the 
versus u dependence with the experimental one. 
On the other hand, as seen in Figs.j^, ^, and 0, all calculations give at high frequencies the 
low-lying Routhian 7r[5, 3, 0]l/2(r=+z) that might produce configurations which are lower 
in energy, cf. Fig.|^, bottom part, than those based on excitations to the tt[6, 6, 0]l/2(r=+i) 
state (cf. also Fig. 0, bottom part). 



2. Bands ^^^Tb{5) and ^^^Tb{6) 

The experimental results for bands ^^^Tb(5) and ^^^Tb(6) indicate a correlation with 
the yrast sequence of the isotope ^^°Tb |2^. This suggests that the neutron Routhian 
z/[7,7,0]l/2(r=— i), cf. Figs. and |^, is not occupied in the related configurations. [Due to the 
level interactions at large cu- values the label of this level changes into z/[7,6,l]3/2(r=— i)]. In- 
deed, calculations indicate that several neutron particle-hole excitations give the required de- 
pendence of JT"*-^-* versus to, as seen in Fig.|^. These are excitations from the z/[7,6,l]3/2(r=— i) 
orbital to either of the i/[4,0,2]5/2(r=±z), i/[5,2,l]3/2(r=±z), or z/[5,l,4]9/2(r=±i) or- 
bitals. A comparison of the dynamical moments corresponding to the configuration 
z/[7,6,l]3/2(r=— z) — > z/[5,2,l]3/2(r=±i) is given in Fig. ^ as an example, while the excita- 
tion energies are presented in Fig.|^. One can see that the energies of the latter configurations 
are ~ (600 to 800) keV higher than those of the former ones, thus favoring the involvement 
of the neutron z/[4,0,2]5/2 and i/[5,2,l]3/2 orbitals as a more likely explanation. 



3. Bands ^^^Tb{7) and ^^^Tb{8) 

The experimental results for bands ^^^Tb(7) and ^^^Tb(8) indicate a correlation with 
the yrast band in the ^"^^Gd nucleus p2|- Such a correlation suggests that the last proton 



intruder level, 7r[6,4,2]5/2(r=— i), remains unoccupied with the corresponding proton on a 
low lying orbital with a small ^T'^^-'-contribution. In fact, the low-lying vr [5,3,0] l/2(r=-|-i) 



17 



orbital visible, e.g., in Fig. |l| provides an excellent result, cf. Fig.|TT]. However, since a pair of 
degenerate bands is seen in the experiment, one should also consider some other possibilities. 
In fact, several neutron configurations mentioned in the preceding Section remain equally 
good candidates, and in Fig. |lT] we also present a comparison with the experiment of the 
results obtained for the i/[7,6,l]3/2(r=— z) z/[4,0,2]5/2(r=±i) excitations. 

It is important to observe that the similarity between the jT'-^-'-moments in ^'^^Gd(l), 
and those in the ^^^Tb(7,8) bands is accompanied by noticeable differences in both the 
Q20 and Q40 moments calculated for the bands with the configurations mentioned. For 
instance the proton excited configuration carries about 0.5 b more than the neutron one: 
Q2ok(651)3/2 ^ 7r(530)l/2] > Q2o[iy{761)3/2 i/(402)5/2], as it is seen in Fig.|lT. 



For the bands ^^^Tb(7) and ^^^Tb(8) we thus have three possible configurations which 
give almost indistinguishable dynamical moments. The importance of the above example 
lies in the fact that the configurations for rather different quadrupole moments and/or de- 
formations may still lead to very similar dynamical moments. Therefore, a direct association 
of "identical" quadrupole moments and "identical" dynamical moments is not supported by 
microscopic calculations. The bands having identical dynamical moments may still carry 
quite different deformations and/or quadrupole moments. 



C. Predictions related to spins in the superdeformed bands of ^^^Tb nucleus 

The problem of the exit spins in SD bands is one of the most actual in the present- 
day studies of the nuclear superdeformation. Many questions remain half-answered since 
the connections between the SD bands and the rest of the nuclear decay schemes are not 
known in the y4~150 nuclei from the experimental point of view. This makes it impossible 
to measure, e.g., the absolute alignments or i7^^^-moments, and to estimate the degree to 
which the rapid nuclear collective rotation resembles that of a rigid body. Such questions 
are being posed also in relation to the pairing behavior at high spins and in many other 
contexts. 

Our calculations allow, in principle, to predict the absolute (and thus also the relative) 
exit spins of the SD bands studied. In practice, however, there are some deficiencies of 
(all) the methods used which impose limitations on the predictive power related to these 
particular results of calculations. 

First of all, in our HF calculations the pairing correlations have been neglected and the 
obtained results concerning the alignment effect may introduce at least small uncertainties. 
Although it is often believed that the spins calculated at a given cranking frequency u must 
be lower when the pairing correlations are included, the actual tendency will depend on the 
position of the Fermi level relative to the occupation of the large-j shells: the down-sloping 
(beginning of the shell) and the up-sloping (upper part of the shell) orbitals contribute to 
the total alignment in opposite ways. Several aspects related to pairing will be presented in 



Sec. VE 



We have analyzed the calculated spin versus frequency plots for all the bands studied in 
this article. The experimental spins of the last levels seen experimentally are not known. 
Therefore, we can only compare the theoretical and experimental differences of spin along 
the calculated and observed bands. To this end, we define the theoretical difference of spin 

as 
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= r'^(cu) - r'^(cuo) + 2h, (5.5) 

where Uo is calculated from the two lowest experimental transitions according to the standard 
expression ( |2.1aD . This theoretical difference can be compared to the experimental difference 



of spin 

^cxpj _ jcxp _ jexp ^2h = 2nh, (5.6) 

(as a function of a;) relative to the (unknown) spin If^^ of the state fed by the lowest 
transition seen in experiment. This experimental difference of spin simply equals twice the 
number n of the gamma transition counted from the bottom of the band. A shift of 2h is 
added to both definitions ( p.5| ) and ( [5.6| ) to account for the fact that the lowest experimental 
value of the frequency (|2.1a|) (ujo) is known at I^^^+2h. 

For all the SD bands in ^^^Tb, calculated within the HF method using the SkM* force 
and with the scaling factor /=0.9 discussed in Sec. |V A 2| , the differences between A^^I and 
^expj j^Q^ exceed Ih, for all the frequencies studied, i.e., from 0.3 to 0.9MeV /h. The 
same is also true (up to 0.5h) for the yrast band in ^^^Dy. On the other hand, similar 
comparison gives a disagreement for the yrast band in ^^°Gd, which indicates a necessity to 
consider the pairing correlations in this nucleus, see Sec. |V h] . 

The spins calculated at Uo give predictions for the exit spins Jo in the SD bands. 



2h, (5.7) 



where the brackets [] denote the rounding to the nearest integer or half-integer compatible 
with the total signature and the parity of the number of particles. The obtained values of 
Jo are collected in Table |I| together with the summary of the theoretical HF configurations 
used. In some cases where clearly more than one configuration may be considered, two of 
them are included in the Table to illustrate the possible uncertainty of the configuration. 

The results for the Jo-values have been given in the Table as differences with respect to 
the J^^-'=69/2 ^ value obtained for Band 1 which is used as the reference. The representation 
in terms of differences is used to stress that the absolute value of J^^-* may be considered less 
certain (given the uncertainties of the scaling procedure, neglect of pairing, etc.) as compared 
to relative alignments between different bands, which are expected to be much less model 
dependent. For example, the exit spin obtained in ^^^Dy(l) from the HF calculations without 
pairing is Io=26h (for the state fed by the 602.4 keV transition). This is fully compatible 
with the exit spin of 57 /2h obtained in band ^^^Tb(2) (for the next transition of 646.5 keV) 
and giving the relative alignment of +1/2 h, cf. Ref. 

The effects of dynamical correlations caused by pairing on the alignment mechanism 
in SD bands of the neighboring nucleus ^^^Dy have been studied, e.g., in Ref. |jl6[. These 



effects are expected to be small, of the order of a few h at the most, and to be rather regular 
functions of the frequency. The corresponding alignment corrections were calculated to be 
negative for several bands in the yl~150 nuclei. We have repeated the analysis of these 
correlations for the even-even neighbors of ^^^Tb, i.e., for ^^^Dy and ^^°Gd, for which the 
pairing effects are not influenced by the blocking mechanism, and are easier to extract in a 
relatively pure form (cf. Sec. |V ii| ). 

Without pairing, the same value of Io=26h is for ^^^Dy(l) obtained both within the 



HF and WS calculations. As shown in Sec. V_E, the inclusion of pairing in the WS method 
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brings this value down to 24 h. Because the pairing influence increases from ^^^Dy towards 
^^''Gd, we can very tentatively estimate that the HF absolute exit spins in Terbium obtained 
without pairing overestimate the values with pairing by about Ah. Hence the inclusion of 
pairing would give the exit spin of I^^^=61/2h for the state fed by the 768.4 keV transition 
in the ^^^Tb(l) band. 

It is important to observe that the obtained exit-spin values vary rather markedly from 
band to band. In particular, in ^^^Tb(l) the value of Jo is signiflcantly larger (by 6^) than 
that predicted for the flrst excited band ^^^Tb(2). 



D. Shape evolution with spin 

The shape evolution for various conflgurations studied in this article has been deter- 
mined by calculating the expectation values of the quadrupole moments, Q20 and Q22, and 
the hexadecapole moments Q40, Q42, and (^44 The calculated axial quadrupole and 



hexadecapole moments, Q20 and decrease with increasing rotational frequency for all 
the bands studied. The corresponding results for the most important proton moments, i.e., 
Q20 and (540 are presented in Fig. |T2|. In Fig. |T3| we show for each band the proton quadru- 
pole moments in function of the dynamical moments. One can see that these relations are 
nearly linear, which indicates that the changes in are correlated with changes of Q2o- 
One can interpret this characteristic result as a dynamical-moment evolution accompanying 
the nuclear-shape polarization due to the rotation. 

It is also worth observing that the proton excited conflgurations have noticeably higher 
quadrupole moments and that their slopes are slightly higher than those of the neutron 
excited bands. 

Slight but systematic deviations from the axiality are predicted with increasing rota- 
tional frequency. It is worth mentioning that the deviations from axiality for the signature- 
partner conflgurations, 7r[3,0,l]l/2(r=+i) 7r[6,6,0]l/2(r=+i) and 7r[3,0,l]l/2(r=— i) — >• 
7r[6,6,0]l/2(r=+i), attributed to Bands 2 and 4, respectively, have opposite signs, although 
they both remain small. 

The axial hexadecapole moments Q4Q are relatively diversifled in that their variations 
in function of the conflguration and rotational frequency may reach as much as 25%. The 
more exotic hexadecapole moments with the non-zero /x-components, i.e., Q42 and Q44, have 
also been calculated. The interest in those moments has grown recently in relation with 
the so-called A/=2 staggering observed in transition energies of some rare-earth SD nuclei. 
The staggering could manifest a presence of those exotic multipoles in the nuclear shapes 
and in particular, the non-zero Q44-moments could manifest the nuclear C4-symmetry, as 
suggested in Ref. ||5^ and discussed in [Q. The effect in question may be very small, as 



shown by microscopic calculations [p^58||. Experimentally, the staggering phenomenon has 
not been reported to exist in ^^^Tb nucleus. Nevertheless, we found it of interest to estimate 
the order of magnitude of the relevant hexadecapole moments. The calculations indicate 
that the Q42-moments, although typically two orders of magnitude smaller than the Q40- 
moments, remain non-zero, while the (544-moments are still at least by a factor of a few 
smaller. 
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E. The moments of inertia in the presence of pairing 



The main purpose of the scahng procedure introduced in Sec. |V was, on the one 



hand, to quantify the systematic deviations from the experiment of the obtained results 
for the jr*^^)-moments. On the other hand this procedure proves indeed that the obtained 
discrepancies are systematic since the same scahng factor has been successfully used in 
various nuclei (and for all the bands known in each of them). The comparisons were made 
for several dozens of bands; presentation of these results was limited so far in this paper 
to 8 bands in ^^^Tb only. Discrepancies of similar order of magnitude were also present in 
calculations of other authors, but were not discussed in any major detail. Only as a result 
of an extended comparison with the experimental data we were lead to a suggestion that 
the discrepancies are systematic and thus that the underlying theories (or parametrizations 
of the forces) might be, and most likely should be improved. 

As it is well known, the pairing correlations influence the moments of inertia in a rather 
characteristic manner. Typically (although not generally), the stronger the pairing correla- 
tions, the smaller the kinematical moments J^^^ . The purpose of the present Section is to 
estimate the effect of pairing on the discussed moments, thus putting more restriction on 
the uncertainties produced by the applied nuclear effective forces. 

A technique used here differs from those of Refs. []T3| and in that we perform the 
particle number projection (PNP) of the rotating BCS type wave functions before variation 
(in the former publication the RPA technique within the Nilsson model has been used while in 
the latter, an approximate projection in terms of the Lipkin-Nogami method together with 
the HFB approximation). Here we apply the pairing-self-consistent Bogolyubov method 
and the PNP before variation, using the standard monopole pairing Hamiltonian with the 
constant average matrix elements, different for the protons and the neutrons (see below), 
and the WS cranking approach. 

Calculations of the SD bands in y4~150 nuclei have shown that the inclusion of pairing 
through the density-dependent zero-range interactions improves the agreement with data in 
^^°Gd, but at the same time the agreement in ^^^Dy is destroyed. It seems therefore, that 
at present one needs some more effort, also in the case of the HF approach, to clarify the 
problem of the effective interactions used, and although here we will not be able to carry 
through such a complex project, we would nevertheless like to quantify the difficulty and 
prepare the starting point for a more advanced discussion. 

Below we are going to concentrate on the two neighbor-isotones of the ^^^Tbse nucleus, 
i.e., ^^^Gdge and ^^^Dygg. We beheve that studying the two even-even nuclei has, in the 
present context, some important advantages. As our calculations, as well as those of other 
authors indicate, in the isotones in question the main pairing contributions come from the 
A^=86 neutron system. Since the relative proton/neutron pairing strengths are uncertain, in 
studying the weak pairing limit, the two even-even nuclei mentioned offer less uncertainties 
as compared to odd-A ones in which the blocking in the case of the odd proton may cause 
an additional difficulty. In other words: it is not known at present whether at the weak 
pairing limit that we are confronted with at high spins the usual blocking procedure is a 
sufficient an approximation. 

In order to estimate the influence of the pairing on the alignment process, we have cal- 
culated the kinematical and dynamical moments in function of rotational frequency at fixed 
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deformations representing the average equilibrium-deformation values for the two nuclei in 



question (for details see captions to Figs. 14 and |l^). We have taken and Z single-particle 
states into account when constructing the Bogolyubov transformation for neutrons and pro- 
tons, respectively. When parametrizing the pairing constants (average matrix elements) we 
took as a starting point the pairing force strength parameters from Ref. ||5^ 

= 1 [18.95 - 0.078(A^ - Z)] , (5.8a) 

Gp = 1 [17.90 + 0.176(iV- Z)] , (5.8b) 

where they have been adjusted to experimental data on normal nuclei by applying the 
pairing-self-consistent Bogolyubov formalism without the PNP. When used within the par- 
ticle projection formalism, these values have to be decreased, by about 15% (several aspects 
related to the PNP technique can be found in Ref. [|^ (solvable model) and, e.g., in Refs. 
I8| , p!6[] (cranking with the WS Hamiltonian)). In the present application we use the reduc- 
tion factors of 0.85 for both the neutrons and the protons, i.e., we apply 0.85 Gn and 0.85 Gp 
in the monopole pairing Hamiltonian. In this respect we in fact apply the parametrization 
studied earlier in the normal nuclei thus not introducing new parameters for those nuclei 
studied in this paper. The results corresponding to our PNP calculations are presented in 
Figs. and |15| for the ^^^Dy and ^^'^Gd nuclei, respectively. Let us observe that the high 
spin results for the dynamical moments in ^^°Gd agree very well with experiment. Note that 
the scaling factor used here, /=0.9, is the same as that applied in the calculations without 
pairing presented previously. In other words, in this nucleus the calculations without pairing 
would not agree with the data. 

Altogether, similarly as in Ref. [0 the results with pairing correspond to an increase in 
j7'(2) moments as compared to those without pairing. At lower spins the crossing observed 
in the data through the rise in the JT*-^-' moments appears at a similar frequency also in the 
calculations, although the theoretical crossing frequency is slightly too large. In any case 
such a crossing is to a far an extent a matter of individual level positions rather than the 
collective pairing phase, and as such is of less importance for the discussion in this Section. 

The results for the dynamical moments in ^^^Dy present similar defect at small frequency 
limit where the calculation provide a crossing-like structure at frequencies which are slightly 
too high. Again, for most of the spin values, the dynamical moments obtained with pairing 
are larger than those without pairing. Consequently, the pairing cannot be a remedy for the 
too high moments as calculated (without scaling) both within the WS and HF methods. 

Concerning the high spin behavior of the dynamical moments (slightly different slopes 
in theory as compared to the data): we have checked that a shape evolution (ignored in the 
calculations presented in Figs. ^ and ^5|) brings in the effects of broad proton crossings whose 
onsets are already visible in the Routhians from Fig.|^. These effects are present in a form 
of a slight modification of the J^'^^ vs. uj dependence in the experimental frequency range 
at highest spins only when the proton pairing is taken into account. The aforementioned 
modification is sufficient to provide the correct asymptotic of the JT^^-* at high spins in ^^^Dy. 

Having obtained an overall agreement with the data on the dynamical moments we have 
also compared the jT'^^^-moments as well. The agreement of the kinematical moments can 
be considered very good for both nuclei when the experimental exit-spin values (so far not 
measured) are assumed to be equal to those calculated according to the method described in 
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Sec |V C| , i.e., /o=24 and 32 h for ^^^Dy and ^™Gd, respectively. Incidentally, without pairing, 
the same method gives the value of /o=38 h in ^^''Gd, and an incorrect behavior for the spin 
difference A^^I. 

Let us note that by taking into account the pairing correlations a good description of the 
kinematical moments at the high-spin limit is obtained for the same values of the scaling 
factor and pairing constants common for both nuclei, although the effects of pairing, espe- 
cially in the neutrons, are calculated to be markedly different in the two isotones discussed. 
In all compared situations a scaling factor was necessary in order to bring the too high no- 
pairing results on the J'^^^ and jT'^^-'-moments towards the experimental data simultaneously 
for both nuclei studied. 

For all the three nuclei, i.e., ^^^Gdse, ^^^Tbse, and in ^^^Dyge, the inertia contribution 
coming from the neutrons is expected to dominate that of the protons by about a factor of 
two; also the pairing effects are calculated to have about a factor of two more effect in the 
neutron case. Comparison of the results in the top and the bottom parts of Figs.|T^ and H 
indicates that at the highest frequency range the differences between the proton contributions 
for the case with pairing as compared to the no-pairing case amount to not more than about 
l^^MeV"^ (about 1%) in ^^^Dy nucleus and about O.bh'^MeV-^ (about 0.5%) in the ^^^Gd 
nucleus. In the case of the neutrons the corresponding differences are slightly larger but also 
very small and thus the pairing is calculated to influence the kinematical moments at the 
highest frequency limit by about 3%. 

At the low frequency range the pairing correlations are expected to be systematically 
more important. The effect of pairing at hu ~ OAMeV on the moment of inertia is about 
10% in ^^^Dy and about 20% in ^^°Gd. The main difference between the calculations for 
^^^Dy and ^^°Gd is the smaller deformation in the case of ^^°Gd; the deformations for most 
of the bands in the ^^^Tb nucleus are in between those for the discussed even-even isotones. 
Consequently we should expect that at the low frequency range in the ^^^Tb nucleus, the 
pairing correlations should still amount to ~ 15% on the average, those of the protons being 
negligible for two reasons: firstly they are very small in the two neighboring even-even nuclei 
and, secondly, the 65th odd proton in the ^^^Tb is expected to diminish the proton pairing 
effect even more due to the blocking arguments. This is also why the results for the Terbium 
nucleus are expected to lie in between those for the neighboring even-even nuclei studied 
here. The detailed analysis of this nucleus calls first, however, for a microscopic justification 
of the scaling factor introduced in the present study, and for an improved treatment of 
pairing in an odd nucleus. 

There exists of course an uncertainty margin in the above estimates. The dynamical 
moments are sensitive to the detailed variation of the single-nucleonic Routhians with the 
rotational frequency and a slow alignment of one single orbital may modify the JT^^^ moment 
considerably thus increasing the risk of estimating the collective pairing effect from the 
influence of one level. Secondly, it is impossible to experimentally separate the proton 
and neutron contributions and thus in the calculations one may obtain in principle similar 
results by increasing, e.g., the proton pairing strength while decreasing that of the neutrons. 
Fortunately, in the two nuclei discussed the protons contribute signiflcantly less to the total 
inertia as compared to the neutrons and, in addition, the neutron pairing effects at = 86 in 
the Dy and Gd isotones differ considerably. Consequently the second uncertainty mentioned 
has a rather small influence on the flnal result. 
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VI. SUMMARY AND CONCLUSIONS 



In this article the two most efficient microscopic approaches in the high-spin nuclear 
structure calculations, i.e., those based on the average nuclear field parametrization in terms 
of the deformed WS potential and the self-consistent HF approach together with the cranking 
method have been applied. The results of the calculations employing these two different 
methods and also the results of the calculations with two different parametrizations of the 
Skyrme interactions have been compared. 

Despite of the fact that the formalisms used differ strongly, the efforts put earher (by 
other authors) in optimizing the parametrizations of the interactions seem to converge: the 
single-nucleon spectra generated by the HF and the WS Hamiltonians are very similar. In 
particular, we have demonstrated that at least for the nucleus studied the self-consistent HF 
solutions and the WS-Strutinsky energy minimized results lead to similar super-deformed 
shell closures. In protons there is a series of shell gaps corresponding to Z=62, 64 and 66 
with a characteristic high level-density areas closing those gaps from below and from above. 
The 7r[3,0,l]l/2 orbital separates the former two gaps and the 7r[6,5,l]3/2 orbital the latter 
two gaps in both SkM* and SHI parametrizations of the Skyrme interaction and in the WS 
case. Of course in the high-level density areas the relative positions of the single-particle 
states differ shghtly but practically all the same Nilsson orbitals are represented there. In 
the neutron spectra the A'"=86 SD gap is cut into two (A'"=86 and N—S5 gaps) by one of the 
signature partners of the z/[7,7,0]l/2 orbital. The neutron spectra manifest two systematic 
differences (when comparing HF and WS): (i) the signature splitting of the z/[7,7,0]l/2 orbital 
is much stronger in the HF approach independently of the force, and (ii) below the big gaps 
we find, in the decreasing order, the z/[6,4,2]5/2, z/[6,5,l]l/2, and i/[4,l,l]l/2 orbitals; in 
the WS case these N=6 orbitals are inverted. The inversion, however, did not cause any 
significant problem in interpretations of the spectra in the ^^^Tb nucleus. 

Theoretical calculations of the dynamical moments indicate a fair or good agreement 
for eight bands known experimentally in ^^^Tb. This is possible, however, only after a 
phenomenological scaling factor, common for all bands, is introduced. Generally, for higher- 
lying SD bands a unique interpretation of the single-nucleonic structure in terms of the 
Nilsson labels does not seem possible (at least at present) for several reasons. Firstly, the 
single-nucleonic Routhian spectra contain levels of different intrinsic characteristics which 
remain nearly degenerate over the long ranges of rotational frequency. They lead to nu- 
merous configurations with very similar alignment and dynamical-moment properties and 
of very close total energies. Secondly, the band degeneracy (twinned band effect) does not 
provide major help in such cases since several among excited bands may produce nearly 
identical band spectra. 

The deformation changes in function of the band-configurations and rotational frequency 
are very characteristic: in all bands discussed not only the quadrupole {Q20) and the hexade- 
capole (Qio) moments decrease with increasing frequency, as observed also by other authors, 
but also their variation seems clearly correlated with that of the dynamical and kinematical 
moments. With progressing alignment small triaxialities (equivalent to the 7-deformations 
of the order of a few degrees) develop as well. 

The possibility of the existence of some special symmetries (like, e.g., an occurrence of 
the shapes with the C4-symmetry) has been studied. No indication for the (544-inoments 
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has been found (an upper limit is predicted at the level of at most 10~^ XQ40). 

The calculations allow to predict the values of the exit spins, i.e., the spins of the lowest 
observed SD transitions. These predictions are given by the values of spin calculated at 
the rotation frequency corresponding to the lowest observed transition, and rounded to the 
nearest half-integer value compatible with the total signature. The differences between the 
exit-spin values in different bands are expected to be more reliable than the corresponding 
absolute values. 

We found that the systematic discrepancies between the experimental and calculated 
results exist in both WS and HF approaches, even if the pairing correlations are taken into 
account. This suggests a need for more effort in determining the optimal parametrizations 
of the effective interactions used for a description of the high-spin properties of the strongly- 
deformed nuclei. 
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TABLES 

TABLE L Experimental 7-ray energies (in keV) in SD Bands 1 and 2 in ^^^Tb. 
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TABLE II. Theoretical configurations of eight bands in ■'■^■'■Tb and the values of exit spins pre- 
dicted within the HF method with the SkM* force and no pairing. The exit-spin values correspond 
to states fed by the last measured 7 transition identified in the last column. The value calculated 
for Band 1 equals to li^^=69/2h. An estimated effect of pairing (see text) would bring this value 
down to 61/2^. Exit spins of all other bands are given with respect to When more than 

one configuration could be a candidate (e.g. in the case of Bands 5 to 7) at most two of them are 
presented in the Table. Spin differences are given in h and the transition energies in keV. The 
second column gives the total parities and signatures of the bands analyzed. 



Band (7r,r) 



Configuration 



2 (-,- 

3 (+,+■ 

4 (-,+■ 

5 (+, - 

6 

6 (-,+^ 



7 
7 
8 



1 (+,-' 




768.4 

-6 646.5 

-3 727.1 

+3 865.8 

-2 710.1 

-2 710.1 

+3 838.2 

+3 838.2 

+1 753.9 

-1 753.9 

785.0 




u[7,6, l]3/2(r 
u[7,6, l]3/2(r 



I 



z/[4, 
i/[4 



0,2]5/2(r 
0,2]5/2(r 



-i) 



30 



FIGURES 



FIG. 1. Single-particle Routhians calculated for the ^^^Tb nucleus by using the HF crank- 
ing formalism with the SkM* interaction. The solutions correspond to the reference configu- 



ration defined in Eq. (5.2). There are four families of curves corresponding to different par- 
ity-signature combinations, namely, the positive-parity Routhians are denoted with the full (r=— i) 
and dash-dotted (r=+z) lines, and the negative-parity Routhians are denoted with the dotted 
(r=+z) and long-dashed (r=— i) lines. Recall that the positive (negative) signature r=+i {r=—i) 
corresponds to the negative (positive) signature exponent a = — ^ (a = +^), the latter notation 
used by some other authors. The curves are labeled using the standard Nilsson labels [A^, n^jAjri 
with the sign of the signature r added at the end. They correspond to dominant components of 
the calculated wave functions at uj= (left set of labels) and a;= 0.8 (right set of labels). 

FIG. 2. Similar to that in Fig.[l|but using the SIII parametrization of the Skyrme interactions. 
The reference configuration has been defined in Eq. (|5.2|). 



FIG. 3. Single-particle Routhians analogous to those shown in Fig.|T| but for the WS cranking 
Hamiltonian. With increasing values of w, the deformation parameters (/92,7,/34) vary so as to 
minimize the total energy (calculated by making use of the Strutinsky method) of the reference 
configuration (|5.2| ). 



FIG. 4. Comparison of the calculated dynamical moments (in Ti^/MeV) for the yrast configu- 
ration in ^^^Tb without (top, /=!) and with (bottom, /=0.9) the scaling procedure described in 
the text. 

FIG. 5. Energy as a function of spin for the proton particle- hole configurations corresponding 
to the low- lying excitations in ^^^Tb. To increase the legibility of the Figure some configurations 
are presented in the top part and other in the bottom part. The labels at the right-hand-side define 
the single-nucleonic particle-hole configurations: the first label gives the hole state, while the one 
following the symbol (— >) corresponds to the particle state. The reference configuration is given 
in Eq. (^.2|). The labels on the left-hand-side give the total parity and the total signature of the 
given configuration. In addition, the configurations are numbered in an arbitrary way, and these 
numbers are used as subscripts to facilitate the identification of configurations. The lines used to 
denote different parity-signature combinations are identical to those used in Fig. |^. The scaling 
has been applied here to conform with the standard of presentation in Fig. ^. 

FIG. 6. Same as in Fig. ^ for the low-lying neutron configurations. 

FIG. 7. The dynamical moments (in ?i^/MeV) of the lowest-energy bands corresponding to 
proton excitations presented in Fig. ^. A presence of at least two types of behavior (two "families" 
of bands) deserves noticing; those presented in the bottom part differ systematically from those 
in the top part. To facilitate the comparison between the curves the vacuum (reference) band has 
been repeated in both, i.e., the top and the bottom parts. 
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FIG. 8. Similar to that in Fig.|^ but for the low- lying neutron configurations. Here at least 
three characteristic forms of behavior are present. The configurations placed in the bottom part of 
the Figure give the slopes which are systematically smaller than those in the top part. Among the 
latter ones the bands demonstrating the crossings in the vicinity of the huj ~ 0.7 MeV are clearly 
visible. 

FIG. 9. A comparison of dynamical moments (in Ti^/MeV) for selected bands presented in 
Figs.0 and ^ with experimental data for Bands 2 (top), 3 (middle), and 4 (bottom). 

FIG. 10. The dynamical moments (in Ti^/MeV) for the neutron configuration proposed for 
Bands 5 (top) and 6 (bottom), compared to the experimental data. 

FIG. 11. Illustration of the possible proton and neutron configurations as those describing 
the structure of Bands 7 (top) and 8 (middle). The dynamical moments (in Tl^/MeV) are nearly 
indistinguishable in proton and neutron configurations, while the corresponding proton quadrupole 
moments (bottom) are markedly different. 

FIG. 12. Evolution of the proton axial moments Q20 and Q40 as functions of the angular 
momentum. The configurations discussed in the text are shown. To gain space the bands are 
labeled by the (vr, r)„-symbols which coincide with those used in previous Figures and are also 
listed in Table ||. 

FIG. 13. Proton quadrupole moments as functions of the dynamical moments (in ?i^/MeV) for 
all eight bands corresponding to these observed experimentally. Note that the high values of the 
quadrupole moments correspond to the low spin limit (cf. also Fig. |l2|) . 

FIG. 14. Comparison with the experimental data of the calculated kinematical (bottom) and 
dynamical (top) moments (in ?i^/MeV) for the yrast configuration in ^^^Dy. Here the cranking 
WS method was applied with deformations (32 = 0.61 and fd^ = 0.12. A small variation of the 
deformation with increasing rotation was not taken into account here since we are interested in 
this case in the effect of the global mechanisms (geometrical scaling, pairing force coupling) similar 
at any deformation. Full lines represent the total quantities, i.e., sums of the proton and neutron 
contributions; labels "t/pair" ("t/no-p") correspond to the pairing correlations being included 
(neglected) respectively. Separate proton and neutron contributions with and without pairing are 
also shown. The scaling factor is /=0.9. The J^^^ moments are obtained from the experimental 
data by applying the calculated exit-spin value of /o=24 for the state fed by the £'^=602.4 keV 
transition. 

FIG. 15. Similar as in Fig.|l4|, but for the yrast configuration in ^^'^Gd. Deformations used: 
P2 = 0.58 and (34 = 0.09. The J^^^ moments are obtained from the experimental data by applying 
the calculated exit-spin value of /o=32 for the state fed by the £'-y=815keV transition. 
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